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Abstract. Basic facts and definitions of conformal moduli of rings and quadrilaterals 
are recalled. Some computational methods are reviewed. For the case of quadrilater- 
als with polygonal sides, some recent results are given. Some numerical experiments 
are presented. This paper is based on [BSVJ and [RVJ. 

1. Introduction 

We give a brief introduction to the conformal moduli of quadrilaterals and rings. 
For a comprehensive survey of this topic see [Kiihj . 

The capacity of condensers has been studied because of its importance in physics 
and its close relation with the potential theory and the theory of conformal and 
quasiconformal mappings. The analytic computation of capacity is possible only for 
very few types of condensers and for this reason several methods have been developed 
for the numerical computation of capacity. 

Let E and F be two disjoint compact sets in the extended complex plane C. 
We assume that each of E and F is the union of a finite number of nondegenerate 
disjoint continua, and that the open set R = C \ (E U F) is connected. Without 
loss of generality, we also assume that oo ^ E. The domain R is a condenser. The 
complementary compact sets E and F are the plates of the condenser. The capacity 
of R is defined by 

(1.1) cap.R = inf / \Vu\ 2 dm, 

u Jr 

where the infimum is taken over all nonnegative, piecewise differentiable functions u 
with compact support in RUE such that u = 1 on E. It is well-known that under the 
assumptions we made above, R is regular for the Dirichlet problem and the harmonic 
function on R with boundary values 1 on £ and on F is the unique function that 
minimizes the integral in (11.11) . This function is called the potential function of the 
condenser. 

Capacity is a conformal invariant: Suppose that / maps R conformally onto R'. 
Let E and F correspond to E' and F' respectively (in the sense of the boundary 
correspondence under conformal mapping). Then capi? = capi?'. This property can 
be used for the analytic computation of capacity provided that the capacity of some 
'canonical' condensers is known and the corresponding conformal mappings can be 
constructed. Unfortunately such an analytic computation can be made only for very 
few doubly-connected condensers; see |ITj . 

If both E and F are connected (and hence R is doubly-connected), R is called 
a ring domain. A ring domain R can be mapped conformally onto the annulus 
{z : 1 < \z\ < e M }, where M = mod-R is the conformal modulus of the ring domain 
R, defined by modi? = 2n/capR. See also [Ahl] . [Hen], jKuhj. 



A Jordan domain D in C with marked points z\,z 2 ,z 3 ,z 4 G dD is a quadri- 
lateral and denoted by (D; z%, Z2, z 3 , Z4) . We use the canonical map onto a rectan- 
gle (D'\ 0, 1, 1 + ih, ih) to define the modulus h of a quadrilateral (D; z±, z 2 , z 3 , Z4) . 
The modulus of (D; z 2 , z 3 , 24, zj.) is We mainly study the situation where the 
boundary of D consists of the polygonal line segments through z%, z 2 , z 3 , z 4 (always 
positively oriented). In this case, the modulus is denoted by QM(£); z\, z 2 , z 3 , 24). If 
the boundary of D consists of straight lines connecting the given boundary points, 
we omit the domain D and denote the quadrilateral and the corresponding modulus 
simply by (z u z 2 , z 3 , z 4 ) and QM(zi, z 2 , z 3 , z 4 ). 





Figure 1: The height of the canonical rectangle is QM(A, B, 0, 1). 

The following problem is known as the Dirichlet-Neumann problem. Let D be a 
region in the complex plane whose boundary dD consists of a finite number of regular 
Jordan curves, so that at every point of the boundary a normal is defined. Let ip to 
be a real- valued continuous function defined on dD. Let dD = A U B where A, B 
both are unions of Jordan arcs. Find a function u satisfying the following conditions: 

(1) u is continuous and differentiable in D. 

(2) u(t) = ip(t), t e A 

(3) If d/dn denotes differentiation in the direction of the exterior normal, then 

^u(t) = ^(t), t G B. 

One can express the modulus of a quadrilateral (D; z\, z 2 , z 3 , z 4 ) in terms of the 
solution of the Dirichlet-Neumann problem as follows. Let = 1, 2, 3, 4 be the arcs 
of dD between (zi,z 2 ) , (Z2,z 3 ) , (z 3 ,z 4 ) , (z&,Zi), respectively. If u is the (unique) 
harmonic solution of the Dirichlet-Neumann problem with boundary values equal to 
on 72, equal to 1 on 74 and with du/dn = on ji U 73 , then by |Ahll p.65/Thm 
4.5]: 

(1.2) QM(£>; z±, z 2 , z 3 , Z4) = / |Vii| 2 dm. 

Jd 

In conclusion, the computation of the modulus of a ring or quadrilateral can be 
reduced to solving the Dirichlet problem (II. ip or the Dirichlet-Neumann problem, 
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respectively. The connection between ring and quadrilateral moduli is given in |Kuhl 
p.102] or [tv] p.36]. 

2. Review of some numerical methods 

Both for the ring and the quadrilateral case we may consider the following 
methods: 

(a) Approximate the canonical conformal map. 

(b) Approximate the solution of the Dirichlet (or Dirichlet-Neumann) problem. 



The recent survey of Wegmann Wcg| provides an extensive review of the var- 



ious techniques of the approximation of conformal mappings. See also Driscoll and 
Trefethen |DrTr] , Papamichael |Pap| . We now mention some of the known methods, 
following closely [BSV] . 

The paper [Gaij of D. Gaier includes a review of the various methods applicable 
to the computation of the capacity of planar ring domains. 

The finite element method was first applied to the computation of capacity by 
G. Opfer [Opf | . Several numerical experiments are reported by J. Weisel | Wei2j . 
Another numerical method is based on the Gauss-Thompson principle which im- 
plies a formula for the capacity involving Green's function. Numerical computations 
are given in |Weilj . N. Papamichael and his collaborators [PKoj . |PS] . |PWa] have 
developed an orthonormalization technique for the approximation of the conformal 
mapping of doubly-connected domains. This technique gives, in particular, approxi- 
mations of capacity. Many numerical computations are presented in the above papers. 

The capacity of a polygonal ring domain can be also computed by the Schwarz- 
Christoffel transformation which provides a semi-explicit formula for the conformal 
mapping of the domain onto an annulus (see [Hen]). For simply-connected domains 
this methods has been developed by T. Driscoll, L.N. Trefethen and their coauthors 
(see |DrTrj . [DrVaj ) . It seems that for doubly-connected polygonal domains the only 
related works are those of H. Daeppen |Dae] and C. Hu [Huj . Hu's method has been 
tested successfully in several computations, (see [Huj . [BV] ). It is partially based on 
the wise choice of certain points on the complementary sets E, F of the ring domain. 

Another numerical-analytic method that can be used for the computation of 
capacity is the multipole method. The potential function is written as a linear 
combination of explicit basic functions (multipoles) with unknown coefficients. The 
coefficients are then computed numerically. The multipoles constitute a complete, 
minimal system in a certain Hardy-type space of functions. Their construction is 
based on the theory of conformal mapping. This method has been developed by 
V.I. Vlasov (see [Via] and references therein) as a general method for numerical so- 
lution of a wide class of boundary value problems. He has applied this approach to 
find the potential function of condensers. 

3. Web-based simulations 

The solutions of the Dirichlet and the Dirichlet-Neumann problems can be ap- 
proximated by the method of finite elements, see [Henl pp. 305-314], [Pap]. Hence, 
this method can also be used to approximate the modulus of quadrilaterals and rings. 
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The Dirichlet-Neumann problem can be numerically solved with AFEM (Adaptive 
FEM) numerical PDE analysis package by Klas Samuelsson. This software applies, 
e.g., to multiply connected polygonal domains. In particular, we may use it to com- 
pute the modulus (capacity) of a bounded ring whose boundary components are 
broken lines. Examples and applications for this software are given in [BSVJ. In 
[HVV] a theoretical formula for computing QM(A, B, 0, 1) was given with its imple- 
mentation with Mathematica. This lead to a more systematic study of the modulus 
of quadrilateral in |DuVuj . In the course of the work on |DuVuj several conjectures 
were formulated and this lead us to look for an improved version of the algorithm in 
[HVVj for the computation of QM(A,B,0, 1). It seems that the AFEM software of 
Samuelsson is very efficient for this purpose. 

Our goal is to use the AFEM software of Samuelsson for computations involving 
moduli of polygonal rings and quadrilaterals, and to write a user-interface providing 
access to AFEM via a web browser (e.g. Mozilla). The advantages in this approach 
are: (1) no programming needed to use AFEM and (2) mobile computing: available 
for everyone. Currently the pilot tests work locally at the local university network. 



Capacity Computation - Mozilla 
File Edit View Go Bookmarks Tools Window Help 



4 , ^ a I! 

Back Forward Reload Stop 



$ http://mat-225.math.helsinki.fi/~rousi/; 



flHome ^Bookmarks -/Red Hat, Inc. ^Red Hat Network CjSupport cSShop H 



Capacity Computation 



** A: Generate inner and outer broken lines by linear representation of a poly gon z |z -> A*z + B*z-con fc 

Outer broken line: 

Number of vertices p in polygon: |5 A: Jl B: [o 

Inner broken line: 



Number of vertices m hi polygon: |3 A: |0. 5 B: |0 

View broken lines only \~ Full summary F submit 

C B: Generate inner and outer broken lines by inserting |x,y)pairs in boxes: 

Outer broken line p-points 
Inner broken line m-points 



FIGURE 2: Entering two regular polygons for computing the capacity of the corre- 
sponding ring domain by using a web browser. 
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Capacity Computation 

1 I , 1 , , , 1 r 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.1 0.6 O.B 1 

Mon flug 08 15:13:27 2005 

Triaii gulLiliun liikI P"[enlial function 



Capacity: 

■ 7.442439e+00 

Figure 3: The output of the program is also displayed by using the web-based 
interface. 

4. Experiments 

We give some examples of problems which can be studied by using the AFEM 
software. 

4.1. Example. Let a,b G C with Ima > 0,lm& > and assume that (a, b, 0,1) 
determines the vertices of a quadrilateral and argfr G (tt/2, it), arg(a — 1) G (0, 7r/2). 
Is it true that 

(4.2) QM(a,fo,0,l) < QM(1 + i\a - 1|, i\b\, 0, 1)? 

1+ila-ll 




1 1 



Figure 4: The quadrilaterals (a, b, 0, 1) and (1 + \a — \b\i, 0, 1) for the left and 
right sides for (14. 2p . respectively. 



Figure 5: The function f(x,y) = QM(1+k, yi, 0, l)-QM(l+a;exp(za), y exp(z/3), 0, 1) 
for x, y G (0,2). Here (3 = 37r/4 and a is n/2 (left) and 7r/8 (right). 



4.3. Remark. The expression on the right hand side of (14. 2p has an analytic expres- 
sion if |a — 1\ = h = \b\ + 1 . Bowman |Bowl pp. 103-104] gives a formula for the 
conformal modulus of the quadrilateral with vertices 1 + hi, (h — 0, and 1 when 
h > 1 as M(h) = X(r)/X(r') where 




= 2h-l. 



Here for < r < 1 

vr X(r') f 1 dx 

Mr) -2XW' %{r) -J T^x^a-r^)' 

and r' = yl — r 2 . 

Therefore, the quadrilateral can be conformally mapped onto the rectangle 1 + 
iM(h), iM(h), 0, 1, with the vertices corresponding to each other. It is clear that 
h-l< M(h) < h . The formula 

M(h) =h + c + 0(e- nh ) } c = -1/2 - log 2/tt , 

is given in [PSj . As fas as we know there is neither an explicit nor asymptotic formula 
for the case when the angle 7r/4 of the trapezoid is equal to a G (0, 7r/2) . 

4.4. Example. (Duplication formula) Let <fi G (0,7r), h, k > 0, A = 1 + hexp(i(j)), 
B = kexp(icj)). We study when the following inequality holds: 

(4.5) QM(A,B,0,1) + QM((1 -B),(l- A), 0,1) < QM(A,B,1 - A,l - B). 

Here equality holds if A = 1 + ih, B = ih and h > . In this special case the result 
may be regarded as a duplication formula. 
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Figure 6: The quadrilaterals in (Oil . 



B 




Figure 7: If B = 1/2 + z(A - 1/2), then A, B,l - A,l - B are the vertices of a 
square and hence QM(A, B, 1 — A, 1 — B) = 1. Letting |A — 1| ->0we see that the 
left side of (14.51) tends to whereas the right side = 1. 



10-, 








Figure 8: The function 

g(x, y) = QM(A B, 1 - A, 1 - B) - QM(A, B, 0, 1) - QM((1 - B), (1 - A), 0, 1) 
where A = x + iy and B = exp(i arg(A — 1)) for x, y G (0, 2). 
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4.6. Remark. Let h, k > 1 and consider next the rectangle with the vertices A = 
1 + i(h + k — 1), B — i(h + k — 1), 0, 1 . If we split this rectangle in two trapezoids 
with the segment joining i(h — 1) and 1 + ih and apply (14. 5p . then we get 

QM(A, B, 0, 1) > QM(A, B, i(h - 1), 1 + ih) + QM(0, 1, 1 + ih, i(h - 1)) . 

If we use the notation from 14. 31 and use the formula QM(A, B, 0, 1) = h + k — 1 then 
we can express this as 

h + k — l> M(h) + M(k) > h + k - 2 . 

4.7. Example. (Open problem |DuVuj ) Fix r,s > 0,a <G (0, 7r/2), /? G (7r/2,7r). 
Determine t > by the condition that the quadrilaterals Q\ — (1 + 2re*°, 2se J/3 , 0, 1) 
and Q 2 = {t + zr, is, —is, t — ir) have equal areas. Is it true that 

(4.8) QM(1 + 2re ia , 2se i/3 , 0, 1) < QM(t + ir, is, -is, t-ir)? 



a 



b 




C 1 




FIGURE 9: Quadrilaterals Q\ and Q2 of Example I4.7I 
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Figure 10: Function 



h(r, s) = QM(t + ir, is, -is, t-ir)- QM(1 + 2re ia , 2se i/3 , 0, 1) 
for a = 7r/4, (3 = 37r/4. 
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